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TECHNICAL MEMORANDUM X- 73361 

NUMERICAL SOLUTIONS OF NAVIER-STOKES EQUATIONS 
FOR THE STRUCTURE OF A TRAILING VORTEX 

SUMMARY 

This report describes the structure and decay of a trailing vortex through 
the numerical solutions of the full, unsteady Navier-Stokes equations. The 
governing equations are recast in terms of circulation, vorticity, and stream 
function as dependent variables, and a second upwind finite difference scheme 
is used to integrate them with prescribed initial and boundary conditions using 
Univac 1108 and IBM 360 computers. The time-dependent numerical solutions 
are carried out within a cylindrical volume that is coaxial with the vortex with 
prescribed boundary conditions on the upstream face, on the axis of symmetry, 
and on the cylindrical outer boundary . 

Different models of the flow are postulated, and solutions are obtained 
describing the development of the flow as integration proceeds in time. A 
p& Ametric study is undertaken with a view to understand the various phenomena 
that may possibly occur in the trailing vortex. 

Using the Hoffman and Joubert law of circulation at the inflow section, 
results obtained with the present method are compared with the experimental 
data of Chigier and Corsiglia on a Convair 990 wing model and a rectangular 
wing. The calculated results compare reasonably well with the experimental 
results on axial and tangential velocity distributions and their decay, axial 
velocity variation along the axis, and core growth. The flow with different 
values of the governing parameters is computed. This parametric study 
greatly increases our understanding of the manner in which the flow dissipates 
downstream. ' - ■ 

Particular difficulty is experienced in getting convergent solutions with 
high values of the swirl parameter . With an exponentially decaying circulation 
distribution and uniform axial flow at the inflow section and at the outer edge of 
the trailing vortex, we get convergent solutions, but the computer time increases 
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excessively as the swirl ratio increases. Also, it is found that the trailing 
vortices persist to long distances downstream before decay. Due to the slow 
development and decay of the vortex core in streamwise direction, vortex 
bursting could not be obtained within the limits of available computer time when 
"uniform" boundary conditions were used. 

To force bursting within a reasonable computer time, an adverse pres- 
sure gradient was imposed on the trailing vortex in the ambient fluid. A burst- 
ing phenomenon was obtained under this condition. A sudden jump in the core 
radius occurred even before the reversal of the flow on the axis. For a given 
Reynolds number, slight changes in the swirl ratio exercised a large influence 
on the bursting bubble. The flow downstream of the bursting rotated like a solid 
body. There is much need to carry out more numerical experiments to under- 
stand better this bursting phenomenon as a function of the governing parameters 
and interpret the sudden increase in the core growth in the light of available 
theories. 


I. INTRODUCTION 


A basic understanding of the various phenomena occurring in a trailing 
vortex is essential to assess the hazards to an encountering aircraft and to 
predict transport and dissipation rates. Such an understanding is also required 
for the safe operation of the aircraft in the terminal areas and for the develop- 
ment of techniques for increasing the dissipation rates. In this investigation, 
an attempt is made to develop an exact procedure of predicting the structure 
of a trailing vortex under various ambient and prescribed conditions. 

An early attempt to solve the problem was made by Newman [1] , who 
linearized the Navier-Stokes equations and neglected the pressure gradient 
term. He obtained analytical expressions for the circulation distribution and 
for the radial and axial velocities. The main contribution of Newman’ s work 
is the analytical derivation of the expressions depicting the exponential decay of 
circulation and axial and radial velocity distributions. Several investigators 
have succeeded in matching experimental data with the tangential velocity profile 
of Newman, using different models of the eddy viscosity. Hall [2, 3] , Kopecky 
and Torrance [4] , and Bossel [5] used these analytical expressions to provide 
initial conditions for the computation of the flow further downstream. Later, 
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Batchelor [6] obtained an analytical solution of the linearized equations with 
the pressure gradient term* He showed that far downstream, the contribution 
from the logarithmic term becomes important. 

Hall [2] simplified the Navier-Stokes equations by quasi-cylindrical 
approximation and integrated numerically the resulting parabolic nonlinear 
equations. At a certain point in the downstream direction, he found large 
gradients in the flow variables which prohibit the computations from proceeding 
further. This indicates the failure of the quasi-cylindrical approximation. 

Hall [3] interpreted this failure as an indication of the onset of a breakdown 
similar to the failure of the boundary layer equations near the separation point. 
This analogy may not be strictly valid because the failure of the quasi- 
cylindrical approximation occurs when the axial velocity is much greater than 
zero. 


Bossel [5] and Mager [7] solved the governing equations with quasi- 
cylindrical approximations by the approximate methods of representing the 
circulation and axial velocity profiles by polynomials. Bossel [ 5] found, in 
his solution, an infinite number of discrete singularities whose corresponding 
critical swirl values depend on velocity and circulation profile shapes. Mager 
[7] developed a momentum integral method for the core of the trailing vortex 
and used third and fourth degree polynomials to represent the axial and the 
swirl velocity distributions, respectively. He also found that for certain values 
of the velocity ratio and swirl ratio, large gradients in the flow cause the 
failure of the quasi-cylindrical approximation. He indicated the onset of 
bursting similar to the one predicted by Hall [3] . He further obtained results 
on either side of the singularity at breakdown. Sarpkaya [8] obtained results 
of vortex breakdown under adverse pressure gradients and correlated his 
experimental data with the theoretical predictions of Mager [ 7] . It may be 
pointed out here that Mager' s [7] results seem to depend critically on the 
degree of the profiles chosen. There is a need to investigate the flow with 
other forms of profiles, particularly the exponential ones. 

Kuhn and Nielson [9] developed basic equations on considerations of 
linear and angular momentum of the vortex and total pressure variation along 
the axis. Using the Hoffman and Joubert [10] law of circulation distribution at 
the inflow section, Kuhn and Nielson computed the flow downstream for gross 
quantities such as the decay of the maximum tangential velocity, the core 
growth, and the axial velocity distribution on the axis. Comparison of their 
theoretical predictions with the experimental data of Chigier and Corsiglia [11] 
seems to be reasonably good. 
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A number of experiments in actual flight, water tank, smoke chamber, 
and wind tunnels have been conducted to understand the structure and decay of 
trailing vortices [11-14], None of these, so far, gives reliable quantitative 
or qualitative information concerning the structure, transport, and dissipation 
of the trailing vortices. 

So far, theoretical work to investigate the structure of the flow in the 
trailing vortex has been done within the framework of a quasi -cylindrical 
approximation of the Navier-Stokes equations. This analysis is essentially 
valid in the core region and ignores mutual interaction between viscous and 
inviscid flow. As such, it can not predict breakdown phenomena with reason- 
able accuracy. 

In this investigation, numerical solutions of the full Navier-Stokes 
equations are used to describe the structure of the flow in a trailing vortex and 
its decay or possible bursting. An unsteady form of the Navier-Stokes equations 
is derived in terms of circulation, vorticity, and stream function as dependent 
variables. With prescribed conditions at the inflow section on the axis and at the 
outer edge of a trailing vortex, the governing equations are integrated by the 
second upwind finite difference numerical scheme on Univac 1108 and IBM 360 
computers. The flow develops smoothly in time from the initial conditions to 
an asymptotic steady state. The flow is essentially laminar. Turbulence is 
introduced in a simplified way through the concept of eddy viscosity. A fully 
developed trailing vortex is simulated by a swirling flow coupled with an axial 
flow. The problem essentially is to prescribe conditions at the inflow section, 
either from experimental data or from some other theoretical considerations, 
and to compute the flow downstream. A parametric study of the structure of 
the flow in the trailing vortex is made and the theoretical predictions are com- 
pared with the available experimental data . 

1 1. GOVERNING EQUATIONS OF MOTION 


A sketch of the coordinate system and the region of integration is shown 
in Figure 1. The Navier-Stokes equations for an unsteady axisymmetric flow 
in nondimensional variables take the following form: 


9T _3_ 

at + az 


(v„r) +-~- (rv r) 

v Z ' r 0r r ' 


i a 2 r a /i ar 
Re az r + r ar \r 8r 
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(v T]) + (v 1 ]) = 
9t 9Z K Z u 9r v r ' 




{f# + ^[r £ (n,) ]} 


9T 

9Z 


( 2 ) 


1 9^6 9 /1 9jfi\ 

r 9Z 2 + 9r \r 9r/ 71 ’ 


(3) 


where r, (p, Z are the radial, azimuthal and axial coordinates, and v , v , 

r (p 

v are the components of velocity in the corresponding directions. 

z 

The variables are made dimensionless as follows: 


U ® Z» r» 

t= — tf V Z = -G ’ rS= C 


v* v’ v’ 

r <£ _ Z 

'r“ C ’ V Vi " V Z U 


_ r» c i/)’ 

r = cv7 ’ 17 = u“ T?t * * = c 2 u 


Here, a prime indicates a dimensional variable. 

Boundary and Initial Conditions 

r = 0 , Z2 Zj : r =r? = ;/) = 0 

r > 0 , Z - Zj : r , t] and i/> to be prescribed (see Conditions 

I and n) . 


5 



J . ; : y : y ; 

r — co , Z > Z t : r-* r 0 , T? — 0 , tp-~ip 0 (4) 

where r 0 and ipo are the asymptotic values of r and ip (see Section IV) 

r ^ 0 , Z = Z 2 : Numerical boundary conditions discussed in 

j Section IV. 

1i '■[ .'i 

Here, Z = Zj and Z = Z 2 represent the inflow and outflow sections of the region 
under consideration. 

At Z = Z j, the following two sets of conditions are imposed: 

Condition I 

(1) Hoffman and Joubert’ s law [9, 10] describes the circulation distribu- 
tion. In our nondimensional variables, it can be written as 



(5) 

( 6 ) 

(7) 


Ki=r 1 V 1 
Kq = r 0 V 0 
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where r 0 , V 0 are the asymptotic values of the radius and the tangential velocity. 


El 

Ki 


= 0.76 - 0.928 In 


(?) 


( 8 ) 


and 


2irK = ~ U C_ |- 
0 4 « L R 


where 


S = wing area 

R = half of vortex pair separation distance at the inflow section 


C = lift coefficient. 
Xj 


From the continuity of vorticity and circulation considerations [9] , 


r = 0. 5625 r, 

1 


(9a) 


r. — • 0. 515 r 
1 k 0 


(9b) 


Here, we need to know rj from experimental data. 
(2) The axial velocity is represented by 


v = 1 - A[1 - 6p 2 + 8p 3 - 3p 4 ] 

z 


( 10 ) 
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where 


p = r/r 2 , r 2 being the distance from the axis at which the axial velocity 
is equal to the ambient velocity (let r 2 = arj where a is a constant). 


A = 


U - - K'r-O 


CO 


U 


= 1 - V, 


CO 


The axial velocity satisfies the conditions 
(a) at p = 0 , v = 1 - A 

where A is positive for a deficit flow and is negative for an excess flow. 


0v, 


Z 

dp 


= o 


3v„ 


a 2 v 


(b) atp = l . -^ = ^ = 0 


Several numerical experiments were made with a third degree profile 
dropping the condition d 2 v / dp 2 = 0 at p = 1. It was found that a fourth degree 

Zi 

profile gave a smoother merging of the flow with the ambient fluid. It is sur- 
mised that higher degree profiles satisfying vanishing gradients of higher order 
at p = 1 may give an even smoother and faster merging. 


( 3) Vorticily rj = 


3v 3v„ 
r Z 


av„ 


az 


3r ~ dr 


= -12A0 2 r(l - 2/3r+ (3 2 r 2 ) 


( 11 ) 


where 


P = 


V2 
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(4) Stream function 


3(?r 2 + — /3 s r 3 - /J 4 r 4 
5 ) 


Condition H 


Here, we consider an exponentially varying circulation distribution as 


follows: 


r = A(l-e ) 


where A and B are constants. In this investigation, we take A = 0. 554, B = 
8. 0 so that the core radius is approximately (0. 4C) . For 


Z — Z j , r>0 , v =A(1-/3Z) 


. 'Ml-,!/.) — 




ri - 0 • 0 


Numerical experiments with various values of A and /3 are carried out. 
We choose A = 0. 75 and j3 = 0.25 so as to force bursting with reasonable com- 
puter time. 
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III. DESCRI PTION OF NUMERICAL METHOD 


A second upwind finite difference scheme was used to integrate equations 
(l) and (2) for r and rj, and the method of successive over-relaxation was used 
for equation ( 3) for ip. This method was used previously by a number of 
investigators [15] , Notable contributions have been made by Torrance [16] , 
Fromm [17] , and Donovan [18] , who investigated the recirculating fluid motion 
driven by the combined effects of a moving wall and natural convection in 
rectangular cavities. Later, Kopecky and Torrance [4] examined the formation 
of eddies in a rotating stream. Roache and Mueller [19] obtained numerical 
solutions for the back step flows with and without splitter plates and for flows 
over square cavities. They discussed the transportive and conservative prop- 
erties of this method and the importance of numerical boundary conditions to 
give accurate and convergent solutions of the Navier-Stokes equations. The 
numerical method used here is essentially the same as that used by Torrance 
together with other investigators in a number of publications. 

The basic method used for equations (l) and ( 2 ) at interior grid points 
is as follows: the radial and axial velocity components are obtained by center 
differencing of the stream function, 
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Here, the subscripts denote the spatial location of a grid point. Circulation 
and vorticity are then advanced in time: 



where the variable f represents either the circulation or the vorticity rj . 

P = 0 I, Q 1 if f is circulation r 

P = 1 , Q = 0 if f is vorticity v . 


Diffusion and swirl ratio terms in equation (20) are represented by the centered 
difference formula as follows: 



and 





2AZ 



It should be noticed that in equation (22) the latest value of circulation 
is used. 
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The advection terms in the second upwind finite difference scheme are 
treated in a different way. They are approximated by the three-point noncentered 
finite differences as follows: 


8 

az 


(v Z f) AZ [ <V z\,j f t, j ' ^z’l-l./i-l.J ] 


£or (i zh,i > 0 md ^z'i-x.j > 0 


(23a) 


or 


AZ [* v Z*i,j £ iU,j " * v Z*i-l,j f i-l, j] 


for (v ) < 0 and (v ) < 0 

**: i>] " 1-1,] 


where 


, <V Z*i+l,i + ^ V Z*i,i 
V Z i,j “ 2 


and 




(v ) + (v ) 

V Vi,i Vi-lii 


(24) 


Similar expressions are used for 8(v f)/ 9r. If (v,_). . and (v_). , . are of 

r Z i,j Z 1-1,3 

opposite signs, a mixed expression consisting of the corresponding terms in 

equation (23a, b) is used for the finite difference approximation. At each 

iteration, the computer checks the sign of (v_) . . and (v ). .at each interior 

Z 1,3 r i,j 

grid point and picks up the corresponding formula for the advection terms. 
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Computations proceed by advancing time (At) . At each time step, new values of 
circulation, F , and vorticity, rj , are obtained in that order according to equa- 
tion (20) . Then, the new stream function 1 p is calculated from equation ( 3 ) by 
the method of successive over-relaxation as follows: 


f: 


(s+l) _ 


(s) 


1,J 


- (i - + 


rn; . + 
i»J 


U) \ 

b 

2 r 

1 . r 1 

(AZ) 2 r + ^ 

/ A \ 2 Ar . v 2 

(Ar) r - — (Ar) 


(AZ)‘ 


(7 S) 


+ ^ 
i-l, ] / 


r + 


Ar 


(Ar)' 


"&) + 

^i, j+1 


,(s) 


r - (Ar) 


(s+l) 

— • ib' 

2 \}-i 


( 25 ) 


Here, ip. . denotes the stream function at the grid point (i, j) after iterations. 

1, J 

The optimum relaxation factor co, = 1. 7 is found adequate for good convergence 

b 

for the given sy stem of grid points. We notice that the most recent values of 17 
and ip are used in calculating the new ip field. At each time step, equation ( 25 ) 
for ip is iterated until the maximum difference at the grid points between suc- 
cessive values of ; p at each time step is < 10 ~ 4 . Velocity components are then 
calculated according to equation ( 18 ) , and the new values of F and 17 are cal- 
culated at the next time step. 


To avoid computational instability, the critical time step at each grid 
point is given by 


crit 


( *Z>U /. . Ar\ 


AZ \ 2 r / 


Ar 


He(AZ)‘ 


1 - 


Ar 

2 r 


Re(Ar)‘ 


-1 


when the mean velocities (v_). . and (v ). . are positive. Similar expressions 

for (At) are derived when (v„) . and (v ) . are negative or they are of 

crit /j \ 9 j T If j 

opposite signs. In actual computations. At = 0.95 x (At) ..is used. For 

cnt 

static stability, no restraints are imposed on spatial mesh sizes. In Figure 2, 
under the prescribed conditions, numerical experiments have been carried out 
for Ar = 0. 015, 0. 025, and a =2.0; and Ar = 0. 035 and a = 2. 5, AZ being equal 
to ( 2 Ar) in each case. Figure 2 shows that under the prescribed conditions as 
the mesh size increases, gross quantities such as core growth, maximum 
tangential velocity decay and axial velocity distribution remain almost the same. 
The slight increase in axial velocity and core growth for Ar = 0. 035 is primarily 
due to increased value of a. Numerical experiments have also been carried out 
with larger numbers of meshes in axial and radial directions. Increasing the 
mesh numbers in r and Z directions does not have appreciable effects on the 
flow structure so long as the prescribed conditions at the inflow section (i. e. , 
a, A and core radius r) are mild and the value of the swirl ratio is low. 

By suitably changing the inflow conditions, the values of the governing 
parameters, and by increasing the number of iterations, it is possible to remove 
ldnks in the azimuthal components of the velocity, oscillations of the axial 
velocity on the axis, and relatively large variations of the flow quantities a few 
grid points inwards from the outflow section. 

As the computations proceed in time, the flow develops smoothly from 
some initial state to the asymptotic steady state. In general, the criterion 
used to determine the steady state is as follows: 


Max 

i*j 


n+1 



n 



10 


-4 


(27) 


It is found that the equation for circulation converges faster than the vorticity 
equation. Also, the Poisson equation for ip takes a lesser number of iterations 
at each time step as the solution converges. 

With an exponentially varying circulation at the inflow section, the num- 
ber of iterations increases excessively to get steady state solutions satisfying 
the previously mentioned criterion. The number of required iterations increases 
further as the swirl ratio is increased. However, the values of the flow quanti- 
ties do not change appreciably after 
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< 0.0005 


(28) 


n+1 n 

Max rj. . -ri. . 

i.i M >iJ 

To save computer time without losing much accuracy, we adopt the pre- 
viously mentioned criterion of convergence when the initial circulation is 
described by the exponential law. 

The basic equations were integrated on the Univac 1108 and IBM 360 
computers of MSFC. Depending on the severity of the parameters and model 
of the flow, typical run times varied from 15 to 125 minutes. 

IV. Dl SCUSSI0N OF NUMERICAL BOUNDARY CONDITIONS 
AT THE OUTFLOW SECTION AND AT THE OUTER EDGE 
OF THE TRAILING VORTEX 


Extensive numerical experiments were carried out with different condi- 
tions at the outflow section and at the outer edge of the trailing vortex to under- 
stand their effect on the detailed structure of the flow. 


At the outflow section, Z = Z 2 , the following conditions were used: 
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or 


(b) v 
r 


l d£ 

r 9Z 


y .ill 

Z r 8r 




8Z 9r 


( 31 ) 


In each case, at Z = Z 2 , derivatives were approximated by a first order back- 
ward finite difference scheme. 


( 2 ) & 
v ' dz 


i.j 


az 



(32) 


v , v,, , and are computed from ip. 

r z 

In general, it is found that the numerical conditions at the outflow section 
do not exercise any significant effect on the rate of convergence and the upstream 
structure of the flow. In the present computations, boundary conditions set in 
(l) (a) are employed. 


Similar numerical experiments have been carried out with different 
numerical conditions at the outer edge of the trailing vortex. Particular atten- 
tion is given to the following conditions on ipi 


( a ) ($) outer e dg e for all i ^ outer edge at the inflow section 


->] = 0*0 


V — OvO, v_ = 1 
r Z 


(b) Aj“ 2 V-i-W 


Here, v^, andt] are approximated from the updated value of ip by a 
first order backward finite difference scheme. 


Condition ( a) simulates the condition of a frictionless moving wall which 
forbids the development of a boundary layer on its surface. Most of the numer- 
ical computations thui> far have been carried out with this condition based on 
the belief that the wall is sufficiently far from the axis and that in a sizable 
region of flow between the axis and the wall a uniform flow exists. As such, the 
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presence of the frictionless moving wall will not affect the development of the 
flow in the trailing vortex. Conditions in (b) update the stream function in each 
iteration, and thus permit the inflow and outflow of fluid at the boundary. This 
simulates more realistically the condition in actual flows. Numerical experi- 
ments were carried out with logarithmic and exponential laws of circulation, 
with boundary conditions at the outer edge prescribed in. (a) or (b) . It is 
found that conditions (b) are highly destabilizing. In several cases, computa- 
tions in successive iterations converged up to a certain extent and then diverged 
or oscillated. Variation of mesh size or their number, changing the time inter- 
val and modifying other conditions at the edge and at the outflow section, did not 
improve the situation. We obtained results that satisfied our convergence 
criterion given in equation (27) in only one case. In Figures 3 and 4, we com- 
pare the results for axial velocity distribution on the axis and the maximum 
tangential velocity variation with distance downstream when conditions (a) or 
(b) are imposed. We find significant differences, both quantitatively and 
qualitatively, in the gross quantities as well as in the detailed structure of the 
flow. If these results are reliable, they indicate that no matter how far the 
outer edge is located from the axis, the prescribed conditions on it exercise 
significant effects on the structure of the trailing vortex. 

With the conditions ( a) at the outer edge, we obtain convergent solutions 
with Hoffman and Joubert* s law of circulation for low values of swirl ratio, and 
with the exponential law of circulation for all values of swirl ratio. In the latter 
case, computation time increased excessively as the swirl ratio was increased. 
Also, the deceleration of the fluid on the axis is slow indicating the need of much 
larger range of integration if we are to obtain bursting through reversal of the 
flow on the axis. Thus, to get more meaningful numerical results, we need a 
larger computer and more computer time. Since our interest is focused on 
tinderstanding the bursting phenomenon, we subsequently imposed an adverse 
pressure gradient in the ambient fluid to accentuate the deceleration of the 
flow on the axis. 

V. DISCUSSION OF RESULTS 


Numerical solutions of full Navier-Stokes equations are obtained for a 
wide range of governing parameters and prescribed conditions. At the inflow 
section, Hoffman and Joubert’ s law of circulation together with other condi- 
tions (as stated in Section n) are prescribed. These conditions at the inflow 
section introduce the following parameters in the problem: 
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1. Radius r* of the core where the azimuthal component of velocity 
reaches its maximum. 

2. Distance from the axis at which the axial velocity matches the ambient 
velocity. 

3. Deficit or excess A in axial velocity on the axis. Other parameters 
appearing in the basic equations are the swirl ratio and the Reynolds number. 

hi the following discussion, we shall try to understand the manner in 
which the flow pattern in the trailing vortex changes as we vary a parameter 
while others remain fixed. 

In Figures 5 and 6 we compute the flow with an excess or deficit in the 
axial velocity on the axis at the inflow section when Re = 2000, swirl ratio = 

0.2, core radius = 0.1, and a = 2. 5. v„(Z , 0) changes from 0.4 to 1. 576. 

Zi 1 

From Figure 5, we find that with an excess axial velocity, the flow on the axis 
decelerates while with a deficit axial velocity the flow accelerates in the down- 
stream direction. With an almost uniform flow at the inflow section, the flow 
remains almost uniform in the downstream direction. This feature of the flow 
is in general agreement with the available experimental data [11-14] . From 
Figure 6, we observe that as the magnitude of the axial velocity decreases, the 
maximum tangential velocity also decreases. There is hardly any change in 
the core growth as a result of the changes in the axial velocity on the axis at the 
inflow section. 

In Figure 7, we compare the variations in the axial velocity on the axis, 

maximum tangential velocity, and core growth as a result of the changes in 

core radius, rj, from 0. 05 to 0. 25 when Re = 2000, swirl ratio = 0. 2, 

v (Z , 0) = 1. 176, and a = 2. 5, We find that as the inflow core radius increases 
Z 1 

from 0. 05 to 0. 25 in steps of 0. 05, the flow accelerates, maximum tangential 

velocity decays slower, and the rate of growth of the core decreases. This 

leads us to believe that a trailing vortex with a larger inflow core radius will 

persist longer. We further observe that the changes in v^Z , 0) and 

Z 1 

[(v ) /V\,] are more pronounced for smaller values of the core radius. 

lv f max 1 

Under the prescribed conditions, there is hardly any change in v^.( Z, 0) when 
the core radius is greater than 0.15. 


In Figure 8, profiles are drawn for the tangential velocity distribution 

as the inflow core radius becomes 0. 1, 0. 2, and 0.25 with Re = 2000, swirl 

ratio =0.2, and a = 2.5. We find that there is no kink in the v -profile for a 

<P 

core radius of 0.1. A kink in -profile appears when the core radius equals 
0.2. The magnitude of the kink increases as the core radius increases to 0.25. 


In Figure 9, we compare v„( Z, 0) and f(v ) /V.,1 when a = 1.0. 

Z cp max 1 

2.0, Re = 2000, swirl ratio = 0.2, inflow core radius = 0. 105, and v (Z , 0) = 

Z 1 

1.176. We find that as a increases, axial velocity increases and 

[(v /Vj ] decreases slightly. An important feature is that the slight 
cp 1 max 

overshoots in [(v ) /VJ in the entrance region with a = 1. 0 smears away 
, „ „ <P max 1 

when a —2.0. 


In Figures 10 and 11, we compare axial velocity and maximum tangential 

velocity distributions with changes in Reynolds number when swirl ratio = 0. 2, 

core radius =0.1, a = 2.5, and v„(Z , 0) = 0.75 or 1.25. Figure 10 indi- 

Z 1 

cates that with an excess inflow axial velocity, v (Z, 0) increases with Reynolds 

Z 

number, while the opposite trend prevails with a deficit inflow axial velocity. 

Larger variations in the entrance region can be corrected by suitably changing 

other parameters. No rational explanation can be obtained for this opposite 

effect of Reynolds number on v (Z, 0) with an excess or deficit inflow axial 

Z 

velocity. Figure 11 presents graphs for [(v ) / VJ with a deficit or an 

f max 1 

excess inflow axial velocity. Reynolds number changes from 5000 to 15 000. 

We find that in all cases the maximum tangential velocity decays slower and 
the core becomes more slender as the Reynolds number increases. The latter 
results show that a trailing vortex with a greater Reynolds number may persist 
longer before dissipation. 

Figures 12 and 13 present graphs for the various flow quantities with an 

excess and deficit inflow axial velocities when Re = 2000, core radius = 0. 1, and 

a = 2. 5. In Figure 12, v (Z , 0) = 1.176 and the swirl ratio takes values of 

Z 1 

0.0, 0.2, and 0.4. We find that the axial velocity decreases and [(v ) /VJ 

(p max 1 

and the core growth remain unaffected as the swirl ratio increases. Figure 13 

indicates that with a deficit inflow axial velocity v ( Z , 0) = 0.75 and swirl 

Z 1 
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ratio = 0.4, v z (Z, o) decreases continuously except for a slight rise in its 

value at the entrance region. For swirl ratio =0.6, v (Z, 0) decreases 

z 

initially and then increases slightly as the outflow end is reached. There is a 

slight decrease in [ (v ) / V.J as the swirl ratio increases from 0.4 to 0. 6. 

cp max 1 

No appreciable change in core growth is noticed with the increase in swirl 
ratio. As the swirl ratio was further increased to 0. 8, we could not obtain a 
convergent solution even after a 70 minute run of the Univac 1108 computer. 
Successive values of 


[At)] 


max 


tj ~V TJ 
new old 



• . I i 

i,j ■■ ■ ' 

oscillated, indicating either convergence was not possible or it would take too 
long. 


This discussion shows that, in general, the swirl ratio has an effect 
similar to the adverse pressure gradient. It can be traced directly in the 
vorticity equation when the swirl ratio is associated with the circulation gradient 
term which is negative so long as the core expands. It plays a crucial role in 
determining whether a breakdown is a possible phenomenon and the manner in 
which a trailing vortex dissipates. Our findings from past investigations are as 
follows: 

1. With Hoffman and Joubert’ s law of circulation, we could not get 
bursting through the stagnation of the fluid on the axis and the consequent 
reversal of the flow. More numerical experiments are probably needed to 
reach this conclusion decisively. 

2. Trailing vortices do not decay like the viscous vortex predicted by 

Newman, viz. (v ) « (i/n/"z) and (r„/ C) « \/~Z. 

<p max 1 

3. Trailing vortices persist too long or too far downstream. 

It may also be noted that the flow downstream depends critically on the 
conditions at the inflow section and the governing parameters. It has a large 
number of parameters which have opposite effects on the flow structure. It 
may be difficult to derive any general conclusions. The previously mentioned 
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parametric study only gives us some idea of how the flow pattern should change 
when one parameter changes and others are kept constant. 

Let us now compare our theoretical predictions with the experimental 
data of Chigier and Corsiglia [12] in 7 by 10 ft tunnel of NASA’ s Ames Research 
Center. These measurements were made with the Convair 990 wing model and 
with a rectangular wing. These data have certain typical features generally 
found in the measurements of other investigators. Distinguishing features of 
the flow on the two wings are as follows: 

1. hi the Convair 990 wing model tests, we observe that a deficit axial 
flow at the inflow section (Z» /C = 6) accelerates to an almost uniform flow at 
the outflow section ( Z’ / C = 12) . Axial profile merges with the ambient flow 
just outside the core. 

2. In the rectangular wing tests, we see that an excess axial velocity on 
the axis at Z’ / C = 4 decelerates to an almost uniform flow at Z* / C = 9. Axial 
velocity attains the ambient velocity at a distance approximately twice the radius 
of the core. 

Chigier and Corsiglia [12] documented the data for transverse and span- 
wise traverses. In certain cases, there are appreciable differences in the 
measurements of the two traverses. These differences are probably due to the 
presence of a secondary vortex besides the primary trailing vortex. The data 
which are markedly affected by the presence of the secondary vortex are not 
included in the comparison. However, some data which bear slight effect of the 
secondary vortex are included in the graphs. At the initial section, an average 
value of the core radius rj, maximum tangential velocity Vj, and the excess or 
deficit axial velocity on the axis is taken from the experimental data for our 
theoretical model of the flow. Kuhn and Nielson [9] found that Hoffman and 
Joubert’ s law of circulation reproduces the experimental data for both wings 
with reasonable accuracy. Using the eddy viscosity concept, they found the 
Reynolds numbers for the flows on the Convair 990 wing model and on the rec- 
tangular wing as 3333 and 2143, respectively. In our theoretical calculations, 
we used the same law of circulation and the Reynolds numbers. For the sake of 
uniformity and to incorporate more points in the core area, we have taken the 
mid-chord C = 1.0. 
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In Figure 14, we compared the theoretical predictions of the core growth, 
the maximum tangential velocity distribution, and the axial velocity on the axis 
with the corresponding experimental data on the Convair 990 wing model. We 
found that the average value of the data is not affected markedly by the secondary 
vortex, and it approximates reasonably well the theoretical predictions. Figure 
15 presents the tangential velocity distributions at various sections in the down- 
stream direction. We observe a clear tendency of smearing the tangential 
velocity profile with the ambient fluid. Also, under the prescribed conditions 
of the governing parameters, there is hardly any change in the tangential 
velocity profile in the outer portion of the flow at the various downstream 
sections. 

In Figures 16, 17, and 18, we compare the present results for the 
overall characteristics of the flow such as the axial velocity distribution on 
the axis, the maximum tangential velocity distribution, and the core growth 
with downstream direction with the corresponding experimental data. Also in 
Figures 19 through 22, we compare the detailed structure of the flow represented 
by the axial and the tangential velocity profiles with the corresponding data in the 
spanwise and normal traverses. Li all cases, the agreement of the theoretical 
predictions with the experimental data seem to be within the experimental error. 

Attempts were made to obtain vortex bursting in the calculations. Within 
Hoffman and Joubert’ s law of circulation, solutions diverged as the swirl ratio 
was increased to 0. 8 with Re = 2000 and an inflow core radius of 0.4 C. We 
subsequently changed the circulation distribution at the inflow section to an 
exponentially varying law, given by equation (13) , while the other conditions 
were kept the same. Convergent solutions resulted, but the number of iterations 
increased excessively. Comparative study indicates that under similar condi- 
tions convergence is faster with the logarithmic law than with the exponential 
law of circulation. With the exponential law, the number of iterations increases 
rapidly as the swirl ratio is increased. 

To save computer time and obtain bursting within the capability of the 
available computers, it was necessary to impose an adverse pressure gradient, 
given by equation ( 14) , at the outer edge of the trailing vortex so that the flow 
decelerates rapidly and bursting through the reversal of the flow may occur. 

In Figure 23, we compare the results of axial velocity distribution on the 
axis of the trailing vortex when the ambient axial velocity is given by v_ - 

A (l - /3Z) . Here, the computations are made with A = 1.0, /3 = 0. 0; A = 1. 0, 
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/?= 0.2; A = 1.0, ( 3 - 0,25; A = 0.75, (3 = 0.25; and A = 0.5, 0 = 0.25 when 
Re = 10. 0 and swirl ratio =0.8. It is quite evident that the greater deceleration 
of the fluid at the outer edge of the vortex leads to larger retardation of the 
fluid on the axis. With a larger deficit axial flow at the inflow section, there 
is larger deceleration of the fluid on the axis. Comparatively larger changes 
take place near the inflow section. For A = 0. 75, /3 = 0.25, and Z = 4. 0, the 
axial velocity at the outflow section is zero; therefore, we impose a strong 
adverse pressure gradient on the trailing vortex. In the other computations, 

A = 0. 75 and /3 = 0. 25 are taken. 

In Figure 24, we compare the axial velocity variation on the axis when 
the Reynolds number changes from 10 to 50 and the swirl ratio is kept constant 
at 0. 8. We find that with increasing Reynolds number, there is larger decelera- 
tion of the fluid. Actually, bursting of the trailing vortex through the stagnation 
of the fluid on the axis and consequent reversal of the flow have taken place. 
Bursting may take place at lower values of the swirl ratio if the Reynolds number 
is increased. This result is similar to the one shown in Figure 10 with the 
Hoffman and Joubert law of circulation. 

Figure 25 presents curves for the axial velocity variation on the axis 

with the distance downstream at swirl ratios of 0. 6, 0. 65, 0. 7, and 0. 8 with 

the Reynolds number equal to 50. At a swirl ratio = 0. 6, v (Z, 0) decreases 

z 

continuously to zero at the outflow section. For other values of swirl ratio, the 
axial velocity decreases faster in the entrance region and then increases to 
attain zero value at the outflow section. This deceleration is greater with larger 
values of the swirl ratio. It is clear that for Re = 50, bursting of the trailing 
vortex has taken place for a value of the swirl ratio lying between 0. 6 and 0. 65. 

Figure 26 shows the growth of the core for the swirl ratios of 0. 6, 0.65, 
0. 70, and 0. 8 and Re = 50. The striking feature of this graph is the nonlinear 
growth of the core up to a certain distance downstream, followed by a sudden 
jump to the outer edge of the region of integration. The core radius remains 
constant thereafter. As the swirl ratio is increased, the core grows faster and 
the jump in the core growth takes place earlier. It is interesting to observe that 
this jump in the core radius takes place even before the bursting of the trailing 
vortex through the stagnation of the fluid on the axis has taken place. This 
abrupt change in the core growth at a particular section downstream shows that 
some sharp change in the flow characteristics is taking place, indicating some 
sort of breakdown in the sense of the Benjamin [20] and Ludwieg [21] theories. 
There is need to investigate the phenomenon of abrupt growth of the core through 
more numerical solutions of the Navier-Stokes equations and comparing them 
with the results of the available theories. 
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Figure 27 shows the variation of the maximum tangential velocity with 

the downstream distance for various values of the swirl ratio and Re = 50. We 

observe that for each value of the swirl ratio the maximum tangential velocity 

decreases with rapidity up to the point where the abrupt change in the core 

growth takes place, and then it remains constant. With the increase in swirl 

ratio, the ratio [ (v ) / VJ decreases faster. Figure 28 presents the pro- 

<p max 1 

files for the tangential velocity at different stations in the downstream direction, 
with a swirl ratio of 0. 8 and Re = 50, As we move downstream, the point of 
maximum tangential velocity moves towards the outer edge of the trailing vortex 
and the profiles become almost linear. This shows that the fluid in the down- 
stream portion rotates like a solid body. 


Figures 29 through 32 present streamlines for swirl ratios of 0.6, 0.65, 
0. 7, and 0. 8 and Re = 50. In Figure 30, we observe the formation of a bubble. 
Figures 31 and 32 indicate that as the swirl ratio increases, the extent of the 
bubble increases in the longitudinal and transverse directions. In Reference 4, 
the authors obtained solutions depicting formation of eddies in a rotating tube. 
The present results differ from them in certain essential features. 


VI. CONCLUDING REMARKS 


Numerical solutions of the full Navier-Stokes equations have been obtained 
for the structure of a trailing vortex. Various phenomena that may occur in the 
trailing vortex are examined by varying the prescribed conditions and the 
governing parameters. Using the Hoffman and Joubert [10] law of circulation 
at the inflow section, the present theory predicts reasonably well the experi- 
mental data of Chigier and Corsiglia [12] on models of the Convair 990 wing and 
a rectangular wing. Use of an exponential law of circulation at the inflow section 
and an adverse pressure gradient at the outer edge of the trailing vortex lead to 
solutions that depict a vortex bursting through the sudden expansion of the core 
and through stagnation and subsequent reversal of the flow on the axis. 

A better understanding of the various bursting processes that lead to the 
dissipation of the trailing vortex is needed. More numerical solutions of the full 
Navier-Stokes equations with different models of the flow should be obtained and 
the results be compared with available theories. 
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